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Abstract 


The magnetic field in the convective zone of the Sun has a random small- 
scale component with the r.m.s. value substantially exceeding the strength 
of a regular large-scale field. For two Majorana neutrino flavors x two helic- 
ities in the presence of a neutrino transition magnetic moment and nonzero 
neutrino mixing we analyze the displacement of the allowed (Am 2 -sin 2 6)- 
parameter region reconciled for the SuperKamiokande(SK) and radiochemi¬ 
cal (GALLEX, SAGE, Homestake) experiments in dependence on the r.m.s. 
magnetic held value b, or more precisely, on a value pb assuming the transi¬ 
tion magnetic moment p = lO^ 11 pb- In contrast to RSFP in regular magnetic 
fields we find an effective production of electron antineutrinos in the Sun even 
for small neutrino mixing through the cascade conversions —► OeR, 

I'eL —* —► v e R in a random magnetic held that would be a signature of the 

Majorana nature of neutrino if v e R will be registered. Basing on the present 
SK bound on electron antineutrinos we have also found an excluded area in 
the same Am 2 , sin 2 28 -plane and revealed a strong sensitivity to the random 
magnetic held correlation length Lq. 
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I. INTRODUCTION 


Recent results of the SuperKamiokande (SK) experiment |]J have already confirmed the 
solar neutrino deficit at the level less Rsk/ Rssm ^ 0.4. 

Moreover, day/night (D/N) and season neutrino flux variations were analized in this 
experiment and in « 4 years one expects to reach enough statistics in order to confirm 
or to refuse these periods predicted in some theoretical models. D/N variations would be 
a signature of the MSW solution to the Solar Neutrino Problem (SNP). The signature of 
the vacuum oscillations is a seasonal variation of neutrino flux, in addition to geometrical 
seasonal variation l/r 2 (f). And signature of the resonant spin-flavor precession (RSFP)- 
solution is the 11 year periodicity of neutrino flux. 

All three elementary particle physics solutions (vacuum, MSW and RSFP) successfully 
describe the results of all four solar-neutrino experiments, because the suppression factors 
of neutrino fluxes are energy dependent and helioseismic data confirm the Standard Solar 
Model (SSM) with high precision at all radial distances of interest (see recent review by 
Berezinsky ||). 

There exists, however, a problem with large regular solar magnetic fields, essential for 
the RSFP scenario. It is commonly accepted that magnetic fields measured at the surface of 
the Sun are weaker than within interior of the convective zone where this field is supposed 
to be generated. The mean field value over the solar disc is about of order 1 G and in the 
solar spots magnetic field strength reaches 1 kG. 

Because sunspots are considered to be produced from magnetic tubes transported to 
the solar surface due to the boyancy, this figure can be considered as a reasonable order- 
of magnitude observational estimate for the mean magnetic field strength in the region of 
magnetic field generation. In the solar magnetohydrodynamics (see e.g. [[|])) one can explain 
such fields in a self-consistent way if these fields are generated by dynamo mechanism at 
the bottom of the convective zone (or, more specific, in the overshoot layer). But its value 
seems to be too low for effective neutrino conversions. 

The mean magnetic field is however followed by a small scale, random magnetic field. 
This random magnetic field is not directly traced by sunspots or other tracers of solar activ¬ 
ity. This field propagates through convective zone and photosphere drastically decreasing in 
the strength value with an increase of the scale. According to the available understanding of 
solar dynamo, the strength of the random magnetic field inside the convective zone is larger 
than the mean field strength. A direct observational estimation of the ratio between this 
strengthes is not available, however the ratio of order 50 - 100 does not seem impossible. At 
least, the ratio between the mean magnetic field strength and the fluctuation at the solar 
surface is estimated as 50 (see e.g. [Q). 

This is the main reason why we consider here an analogous to the RSFP scenario, an 
aperiodic spin-flavour conversion (ASFC), based on the presence of random magnetic fields 
in the solar convective zone. It turns out that the ASFC is an additional probable way to 
describe the solar neutrino deficit in different energy regions, especially if current and future 
experiments will detect electron antineutrinos from the Sun. The termin “aperiodic” simply 
reflects the exponential behaviour of conversion probabilities in noisy media (cf. [f>] , ||). 

As well as for the RSFP mechanism J7 : ] all arguments for and against the ASFC mecha¬ 
nism with random magnetic fields remain the same ones that have been recently summarized 


2 


and commented by Akhmedov (see || and references therein). 

But contrary to the case of regular magnetic fields we find out that one of the signatures 
of random magnetic held in the Sun is the prediction of a wider allowed region for neutrino 
mixing angle in the presence of v e from the Sun including the case of small mixing angles 
(see below section IV). 

Notice that if electron antineutrinos from the Sun were detected at the Earth this would 
lead to conclusion that neutrinos are Majorana particles and this is a very attractive point 
for applying RSFP in regular magnetic held or ASFC in random helds in the Sun. 

The SK experiment provides a stringent bound on the presence of solar electron antineu¬ 
trinos at least for high energy region A phenomenological and numerical analysis of 
v e L —> h e _R-conversion in the solar twisting magnetic holds JT 0 | - |l!| shows that it is possible 
to obtain a noticeable yield of v e R |ll[ consistent with the limit p]. Twisting magnetic 


helds themselves are, however, very specihc and it is hard to explain their existence and 
origin. Therefore, more realistic models of magnetic held in the Sun are necessary and for a 
start we treat here both SNP solution and z/ e -pro duct ion applying, as we hope, a bit more 
realistic model of random magnetic helds in the convective zone of the Sun. 

The random magnetic held is considered to be maximal somewhere at the bottom of 
convective zone and decaying to the solar surface. To take into account a possibility, that 
the solar dynamo action is possible also just below the bottom of the convective zone (see 
0), we accept, rather arbitrary, that it is distributed at the radial range O.7P 0 —l.OP 0 , 
i.e. it has the same thickness as the convective zone. Of course, our model of the random 
magnetic held (see in details below) is a crude simplification of the real situation, however, 
the results are presumably more or less robust. 

For each realization of random magnetic helds we find a solution of the Cauchi problem 
in the form of a set of wave functions u a (f) —\ v a (t) | exp {ia a (t)) obeying the unitarity 
condition for the probabilities P aa if) = 


Peeif) + Pee{t ) + Pfj,fj,(t) + Pp,p(t) -1 j (1) 

where the subscript a equals to a = e for u e L, a = e for zVr, a = n for and a = ft for 
correspondingly. 

In section II the experimental observables are defined through the neutrino conversion 
probabilities P aa (E) at the detectors (r = R,Earth ) accounting for neutrino propagation both 
in the Sun and on the way from the Sun to the Earth. 

In subsection II.A we give general formulas for the neutrino event spectrum in real time 
experiments with measurement of recoil electron energy in the z/e-scattering. In contrast 
to the well-known case of regular helds, we need to find the probabilities P aa (r = t) as 
the functions of a local position in the transversal plane too, P aa (r) = P aa (r,£,r]), because 
longitudinal prohles of the random helds (along r) are generally different in different points 
in the transversal plane even for an instant when parallel neutrino huxes directed to the 
Earth (called here ’’rays”) cross the convective zone (see Fig. |T). 

As a result of different realizations of random helds along different rays the probabil¬ 
ities Paa( r ,{L, r l) are random functions built on randomness of magnetic helds in all three 
directions. 

The same probabilities for left-handed electron neutrinos, P ee (r, £, 77 ), are used in sub¬ 
section II.B to define neutrino hux measured in radiochemical experiments (Homestake and 
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GALLEX + SAGE, in SNU units). 

In main section III we give some physical arguments supporting the random magnetic 
held model implemented into the master equation Eq. (|TTD (subsection III.B.l). After that 
we describe the mathematical model of random magnetic fields (subsection III.B.2). 

Then in subsection III.C we analyze asymptotic solution of the master equation in the 
case of small neutrino mixing. We briefly discuss there many possible analytic issues: in 
particular, magnetic held correlations of finite radius (subsection III.C.2), linked cluster 
expansion and higher moments of survival probability (subsection III.C.3). We give also 
an interpretation of the random magnetic held influence as a random walk over a circle 
(subsection III.C.4). 

In subsection III.D we describe the algorithm of our numerical approach. The main goal 
is to calculate the mean arithmetic probabilities as functions of mixing parameters sin 2 2 6 , 
5 = Am 2 /4E, 


(Paa{&, sin 2 26)) = J f d^driP aa {5, sin 2 26>, £, 77 )/ J j d£dr) , 


( 2 ) 


under the assumption = $f^(£, rj)Xi(E), where A i(E) is the correspond¬ 

ing normalized differential hux, j'^ rnax( - l > \ % (E)dE = 1 and $f ) ' ) (£, 77) is the integral hux 
of neutrinos of kind ”i” (i = pp , Be ,pep , B) assumed to be constant and uniform, 
< h|°' ) (^,??) = = const , at a given distance r = t from the center of the Sun ||T3| . This 

simplifies the averaging in the transversal plane since 


’ £. 2 + t i 2 <Rco 


d£dr7$S 0) (£,77)P aa (£,77) x [J 


ew<R 2 0 


d^dr,}- 1 = # x (P aa ) . 


( 3 ) 


It is worth to note that the above averaging procedure merely rehects the physical prop¬ 
erties of neutrino detectors which response to the incoherent sum of partial neutrino huxes 
incoming from all visible parts of the Sun. At the same time it is well known that the 
most studied statistical characteristics in statistical physics and especially in the theory of 
disordered media are those that are additive functions of dimensions (cf. [14]). Their main 
distinctive feature is that, being devided by corresponding volume, they become certain in 
macroscopic limit, i.e. self-averaged. From definition Eq. (^) it follows that the integral in 
the numerator in the r.h.s. is indeed an additive function of the area of the convective zone 
layer (see also Fig. |l]) and therefore for increasing area one can expect that (P aa ) should be 
self-averaged. Our results below confirm this property ( subsection III.C.3 and section IV). 

In section IV we analyze how these probabilities depend on the r.m.s. held value yj ( b 2 ) 
and argue why electron antineutrinos are not seen in the SK experiment reconciling available 
experimental data for four solar neutrino experiments. 

I 11 final section V we discuss our results comparing two mechanisms, RSFP and ASFC, 
for the same strength of regular and r.m.s. fields and p = 10 ~ ll ps- 

II. EXPERIMENTAL OBSERVABLE VALUES 


If neutrinos have the transition magnetic moment, then in real time experiments like 
^e-scattering in the SK case [jl] the recoil electron spectrum depends on all four neutrino 
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conversion probabilities, P aa (see subsection II.A), while in the case of radiochemical ex¬ 


periments (Homestake, Gallcx, SAGE [ 161) the number of events measured in SNU units 
depends on charge current contribution with left-handed electron neutrinos only, i.e. on the 
survival probability P ee that itself is a function of the magnetic field parameter g,B±, as well 
as of the fundamental mixing parameters Am 2 , sin 2 26 (subsection II.B). 


A. Spectrum and number of neutrino events in ne-scattering experiments 

In a real time experiment one measures the integral spectrum (the number of neutrino 
events per day) 


iV„= E [ Ttt 'dT dN - m 


dT 


(4) 


where experimentalists devide the whole allowed recoil electron energy interval T > T t h into 
bins AW = T k+1 — T k with T kmin = T th . The energy spectrum of events, dN u /dT , has the 
form 

r r d^di]}- 1 x f n _ d£d7iJ2®i 0 \^v) x 


dT 


= N P 


>ew<Ri 

rEraax (^) 


U 2 W<R 2 co 

dcj {M) 

dEe(E)X l (E)(^^(E,m,rj)) , 


(5) 


where N e is the total number of electrons in the fiducial volume of the detector; E min (T) = 
(T/2)(l + yjl + 2 m e /T) is the minimal neutrino energy obtained from the kinematical in¬ 
equality T < T max = 2E 2 /(2m e + 2 E) and E max {i) is given in ]13[. For simplicity the 
efficiency e(E) above the detector threshold T th is substituted by e(E) = 1. 

Notice that we have assumed the same core radii for all neutrino sources (kinds ”i”), 
H = /?b) 

1l core — 1l core ■ 

The averaged differential cross-section {da^lk{E, T\ l£> v)/dT) in Eq. (H), 


da {M \ r d™- 

(^A( E ,m, V )) = 


dT 


n -1 


L J 0 


X 


exp(—(T - T f/2A(T ) 2 )dT 


x 


62 /o a 

, )Ud - 

l jo dl 


exp(—(T - T Y/2A{T f)dT 


( 6 ) 


is given by the energy resolution AT = A (T) in the Gaussian distribution where T' is the 
true recoil electron energy, T is the measured energy, and by the ue scattering cross-section 


for four active Majorana neutrino components |15 

[Ey 

dT 


AAaama = { ,,) [A+ 4(1 - 

+Pe- e (E, f, ri) \g\ + g 2 eL (l - 


+Pp,p,(E, £, 7]) g R + g 2 L (l —) 


T\ 

) 2 

m e T' 

E' 

1 
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T\ 

) 2 
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1 
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Here g e L = £ + 0.5, g /t L = £ — 0.5, gu = £ = sin 2 2 6w — 0.23 are the coupling constants 
in the Standard Model(SM). 

Let us stress that a specific dependence of probabilities on the transversal coordinates 
£ , r) vanishes for a homogeneous regular magnetic field which is a function of the longitudinal 
coordinate r = t only. Therefore for scenarios with regular fields (for instance, in [j7|) the 
spectrum Eq. (|5|) takes the usual form without dependence on transversal coordinates. 

If we assume uniform partial neutrino fluxes, the integration over the transversal co¬ 
ordinates £, rj in the integrand of the event number Eq. (|J) leads to the averaging of 
probabilities, see Eq. (j3|) and subsection III.D. 


The energy resolution in the SK experiment is given by A (T) = A 10 y/T/10 MeV and 
for the Kamiokande experiment is estimated as AT/T ~ 0.2 at the energy T = 10 MeV, or 
A 10 = 2 MeV. This irreducible systematical error becomes even worse near the threshold 
(T t h = 6.5 MeV at the present time and one plans to reach T t h = 5 MeV in 1998). 

In BOREXINO (starts in 1999), where one has a liquid scintillator, it is expected to 
observe approximately 300 photoelectrons ( phe ) per MeV of deposited recoil electron energy. 
This gives an estimate of the energy resolution of the order m 


AT 

-( BOREXINO ) ~ 



0.058 

y/T/MeV ’ 


( 8 ) 


corresponding to 12 % for the threshold ( T t h — 0.25 MeV) and 7 % for the maximum recoil 
energy for Be neutrinos, T max ~ 0.663 MeV. 

In HELLAZ the Multi-Wire-Chamber (MWC) counts the secondary electrons produced 
by the initial one in helium. It is expected to count 2500 electrons at the threshold energy 
T t h — 0.1 MeV , so in that case |L8| 


AT 


(HELLAZ) ~ 


0.02 


JfJTh 


(9) 


or the energy resolution of the order 1.5 % for the maximum for pp neutrinos T max ~ 
0.26 MeV. 

Other experimental uncertainties must be incorporated into the differential spectra Eq. 
(||). We have neglected an unknown statistical error, which we expect to be less than the 
systematical one after enough exposition time, as it will decrease as | ±AN V /N V |~ H 1 ^ 2 . 
We have also neglected inner and external background contributions to the systematical 
error considering them as the specikcal ones for each experiment. 


B. Radiochemical experiments 

The number of neutrino events in GALLEX (SAGE) and Homestake experiments mea¬ 
sured in SNU (1 SNU = 10 -36 captures per atom/per sec) has the form 

r r^max (0 

1 (T 36 ^ a = E * / K(E)aca, ci(E) < P„(E) > dE , 10 

i JE th (Ga,Cl) 
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where the thresholds E t h(Ga , Cl) for GALLEX (SAGE) and Homestake are 0.233 MeV and 
0.8 MeV correspondingly and the capture cross sections aa a , ci(B) f° r gallium and clorine 
detectors are tabulated in JT3[, see Table 8.4. 

For P ee = 1 the SSM predictions for radiochemical experiments are listed in the Table 8.2. 

T3| . For updated theoretical fluxes in BP95 model |L9] we find the mean SSM predictions 137 
SNU for GALLEX (SAGE) and 9,3 SNU for the Homestake experiment. The experimental 
data JR]] pronounce the electron neutrino deficit, 69.7 ± 6 . 7^5 SNU for the GALLEX, 
69 ± lOt? SNU for the SAGE and 2.55 ± 0.17 ± 0.18 SNU for the Homestake experiments. 

Notice that the ratio of the experimental data to SSM predictions above does not depend 
on theoretical uncertainties of the integral neutrino fluxes, in Eq. m ■ Neutrino 

deficit expressed through these ratios @ is 0.509 ±0.059, 0.504 ±0.085 and 0.274 ±0.027 for 
GALLEX, SAGE and Homestake correspondingly where the experimental lcr-errors come 


from yjc7 2 stat ± a\ 
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III. MASTER EQUATION AND ITS SOLUTION 


A. Master equation 


We consider conversions v e L.n —■> v a L,R, a = p or r, for two neutrino flavors obeying the 
master evolution equation 


( VeL ^ 

( 

VeR 


Vy,L 


\ "vR ) 

V 

cos 20, 

S 2 = 


c 2 5 


0 

s 2 <5 


0 

-V e - c 2 S 

s 2 5 


s 2 5 

-pH_{t) 
Vfj, ± c 2 5 
0 


pH+(t) 

s 2 5 

0 


\ 

( ueL \ 


VeR 

) 

\ Vy,R J 


( 11 ) 


the neutrino active-active transition magnetic moment; H±(t ) = B 0 ±(t) + b±(t), H± = H x ± 
iH y , are the magnetic held components consisting of regular ( B 0 ±) and random ( b ±) parts 
which are perpendicular to the neutrino trajectory in the Sun; V e (t) = GF\/2(p(t) / m p ) (Y e — 
Y n /2) and V^(t) = G j pv / 2(p(t)/m p )(— Y n /2) are the neutrino vector potentials for u eL and 
in the Sun given by the abundances of the electron ( Y e = m p N e (t) / p(t)) and neutron (Y n = 
m p N n (t)/p{t)) components and by the SSM density prohle pit) = 250 g/cm? exp(—10.54t) 

m- 

Notice that we did not assume here a twist held model with analogous constructions 


B± = B x ± iB y in off-diagonal entries of the Hamiltonian in Eq. ([Tip These expressions 
are derived from the initial Majorana equations in the mass-eigenstate representation where 
transversal part of the spin-hip term pa ■ B leads automatically to the ” twist-form” in 
the Schrodinger equation above written in havor representation (see |2(J). We have not 
performed the phase transform of the Hamiltonian in Eq. ([ll|) since, in contrast to 


10 , the 


phase <£>(£) in H± = i/^(t)e ±l$ h' ) is the random one that gives uncertainty in the additional 
term $ for the resonance condition [|l0j. Instead of that we have solved Eq. (|LL|) directly 
using computer simulation (see below). 
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B. Solar magnetic fields 


1. Random magnetic fields 


The r.m.s. random component b is assumed to be stronger then the regular one, B, 
and maybe even much stronger than B 0 , (b 2 ) Bq, provided the large magnetic Reynolds 
number R m leads to the effective dynamo enhancement of small-scale (random) magnetic 
fields. 

Let us give simple estimates of the magnetic Reynolds number R rn = lv/v m in the 
convective zone for fully ionized hydrogen plasma (T ^ Ih ~ 13.6 eV ~ 10 5 K). Here 
l ~ 10 8 cm is the size of eddy (of the order of magnitude of a granule size) with the turbulent 
velocity inside of it v ~ 10 5 cm/s. Provided an equipartition between the turbulent kinetic 
energy and the mean magnetic field is suggested, we obtain v ~ va = B /\/ 47 rp, where va is 
the Alfven velocity for MHD plasma, B is the mean (large-scale) magnetic field in convective 
zone and p is the matter density The magnetic diffusion coefficient, or magnetic viscosity, 
u m = c 2 / Air a con d entering the diffusion term of the Faradey induction equation, 


<9H 

~dt 


= rot[v x H] + z/ m AH , 


( 12 ) 


where H = B + b is the total magnetic field (mean field plus fluctuations), is given by 
the conductivity of the hydrogen plasma a con d = ^pi/A-KVep. Here c is the light velocity; 

c jjpi = \J A r Ke 2 n e /m e = 5.65 x 10 4 1 /n)T s" 1 is the plasma, or Lengmuir frequency; u ep = 
50n e /T 3 / 2 s " 1 is the electron-proton collision frequency; and the electron density n e = n p 
and the temperature T are measured in cm -3 and Kelvin degrees respectively. 

ft follows that the magnetic diffusion coefficient u m ~ 10 13 (T/1 K)~ 3//2 cm 2 s~ l does not 
depend on the charge density n e and it is very small for hot plasma T > 10 5 K. From 
comparison of the first and second terms in the r.h.s. of the Faradey equation Eq. (|T 2 D 
we find that v/l n m /l 2 , or u m -C vl ~ 10 13 cm 2 s~ 1 since T/l K 1. This means that 
magnetic field in the Sun is frozen-in. 


Iv/cmfs 1 )(T/1 A') 3 / 2 . A standard 


As the result we obtain R. rn = Ivuj/J( c 2 v ep ) ~ (10 ^ 13 
estimation for R m in the solar convective zone is R m = 10 8 (see ||21|| ). 

The estimation of b/B for the solar convective zone (and other cosmic dynamos) is the 
matter of current scientific discussions. The most concervative estimate, simply based on 
equipartition concept, is b = const B. According to direct observations of galactic magnetic 
field presumably driven by a dynamo, b ss 1.8 B ( [Q). A more developed theory of equipar¬ 
tition gives, say, b = constB In R m (see ^TJ). For our consideration taking this constant as 
4:/r would be more than sufficient. 

Notice that this estimate is considered now as very concervative. Basing on more detailed 
theories of MHD turbulence estimates like b ~ \/R m B are discussed |23|] , [|4|] . Such estimates 


give a free room for a very large values of b. Let us stress, however, that the larger estimate 
of b is accepted, the more difficult is to refer for a dynamo as an origin for B, so the estimates 
of Vainstein and Cataneo, [p3l and [0 are hardly comparable with the dynamo nature of 


large-scale solar magnetic field. 

Being interested in neutrino conversions, we need only a partial information concerning 
the small-scale magnetic fields. In particular, because of the very rapid propagation of 

















neutrino in comparison with MHD timescales, we are interesting only in a b distribution at 
a given instant in a given direction. 


2. The model of the magnetic field 


The assumed profile of a random magnetic field is presented in the Fig. [I], 

We also show there the profile of the regular magnetic field used in calculations. We find 
that RSFP results are not sensitive to the profile of a regular held. However, switching on 
random fields, we find out an essential change of the probabilities P aa {t) fc> r the same mixing 
parameters sin 2 20, Am 2 /2E. This difference of the magnetic held model influence on the 
solution of the master equation Eq. (Tl) becomes the more significant the larger the r.m.s. 
value b is assumed. 

The numerical implementation for the random magnetic held has been chosen as follows. 
We choose the correlation length for the random held component as L 0 = 10 4 km that is 
close to the mesogranule size. 

Then we suppose that all the volume of the convective zone is covered by a net of 
rectangular domains where the random magnetic held strength vector is constant. The 
magnetic held strength changes smoothly at the boundaries between the neighbour domains 
obeying the Maxwell equations. Since one can not expect the strong inhuence of small details 
in the random magnetic held within and near thin boundary layers between the domains 
this oversimplified model looks applicable. 

In agreement with the SSM |L3| we suppose that the neutrino source inside the Sun 
is located inside the core with the radius of the order of R u = 0.1 R@ = 7 x 10 4 km and 
neglect for the sake of simplicity the spatial distribution of neutrino emissivity from a unit 
of a solid angle of the core image. However, different parallel trajectories directed to the 
Earth cross different magnetic domains because the domain size Lq (in the plane which is 
perpendicular to neutrino trajectories) is much less than the transversal size = R v of the 
full set of parallel rays, L 0 -C R u . The whole number of trajectories (rays) with statistically 
independent magnetic fields is about R 2 /Lq ~ 50. 

Thus, the perpendicular components b x ^ y (r) along one representive ray are the random 
functions with the correlation length Lq (see Fig. ||). 

At the stage of the numerical simulation of the random magnetic held we generate the 
set of the random numbers with the given r.m.s. yj (b 2 ) and suppose that the strength of 
the random held is constant inside the rectangular volume with the radial size Lq and with 
the same diameter L 0 in the £ , 77 -plane. We assume that yj (bf) = J ( b 2 ). 

We generate random values of b in different cells both along the neutrino trajectory and 
in the transversal plane £ , 77 , solve the Cauchy problem for 50 rays and calculate the mean 
probabilities Eq. (||). 

Notice that we substitute into the integrand of the neutrino event number dN v /dT 
the probabilities P aa ( r = Reo,£,,v) a t the Earth taking into account also vacuum neu¬ 
trino oscillation probabilities (for zero magnetic helds and zero density in the solar wind) 
and neglecting the Earth effects for relatively small Am 2 ~ 10 ~ 8 eV 2 ( l v = 2.48 x 
lO~ 3 E(MeV)/Am 2 (eV 2 ) km l 0 = 3.28 x 10 4 km/p(g/cm 3 ) 

for p ~ 4.5-11 g/cm 3 in the Earth). As for MSW-region Am 2 ~ 10 -5 eV 2 we do not 
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consider here the possible Earth influence which could result in day-night variations. 


C. Asymptotic solution of the master equation 


Let us consider how random magnetic field influences the small-mixing MSW solution to 
SNP with the help of some simplified analytical solutions of the master equation Eq. ( 0 . 
Indeed, it turns out that for the SSM exponential density profile, typical borone neutrino 
energies E ~ 7 -j- 14 MeV, and Am 2 ~ 10 -5 eV 2 , the MSW resonance (i.e. the point 
where V e — = 2c 2 5) occurs well below the bottom of the convective zone. Thus we can 

divide the neutrino propagation problem and consider two successive stages. First, after 
generation in the middle of the Sun, neutrinos propagate in the absence of any magnetic 
fields, undergo the non-adiabatic (non-complete) MSW conversion u eL —> is^ L and acquire 
certain nonzero values P ee end P MM , which can be treated as initial conditions at the bottom 
of the convective zone. For small neutrino mixing S 2 —■► 0 the (4 x 4)-master equation 
Eq.( 8 ) then splits into two pairs of independent equations describing correspondingly the 
spin-flavor dynamics v e i —> and —■► v e R in noisy magnetic fields. In addition, once 

the MSW resonance point is far away from the convective zone, one can also omit V e and 
Vn in comparison with C 28 . For u e i —> conversion this results into a two-component 

Schrodinger equation 


( v eL \ ( -8 nb + (t ) 

V/^-W 5 


Al \ 


(13) 


with initial conditions |zz e L| 2 (0) = P ee ( 0 ), |h /t ^| 2 = P^( 0 ) = 0 . As normalized probabilities 
P ee (t ) and Pjifiit) (satisfying the conservation law P ee {t ) + P^it) = P ee (0)) are the only 
observables, it is convenient to recast the Eq.([13]) into an equivalent integral form 

t t\ 

S{t) = 5(0) - 4 1 ? J dp J dt 2 S(t 2 ) 

0 0 

{[b x (ti)b x (t 2 ) + b y (ti)b y (t 2 )} cos25(ti - t 2 ) 

-[bx(ti)b y (t 2 ) - b x (t 2 )b y (ti)] sin25(fi - f 2 )}, (14) 

where S(t) = 2 P ee (t) — P ee ( 0 ) = P ee (t) — Pp,p.(t ) is the third component of the polarization 
vector S. 


1. 8-correlations 

If Lq is the minimal physical scale in the problem, we can consider random magnetic 
fields as ^-correlated: 

2 

< b x {p)b x {t 2 ) + fry(P)fry(£ 2 ) > = - < b 2 > -L 0 5{ti - t 2 ), 

< b X by > = 0. 

In this case it turns out that averaging of the Eq. (|T4D over random magnetic fields is exactly 
equivalent to averaging of its solution (cf. ||). The result reads: 
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(15) 


< S(t) >= 5(0)e _rt , 


where < S(t) > is the mean value, and decay factor 

r = < b 2 > L 0 . 


(16) 


If time is measured in units of Lq (recall that c = 1), supposing that neutrino traversed N 
correlation cells, i.e. t = Lq x N, we have 


< S(N) >= 5(0) exp 


—fJ - 2 <b 2 > L j 


%xN 


(17) 


so that N plays a role of an extensive variable like the volume V in statistical mechanics. 


2. Correlations of finite radius 


Let us consider Eq.(|T4|) in a spirit of our numerical method. Dividing the interval of 
integration into a set of equal intervals of correlation length Lq, for a current cell n we have 

S n = S(n ■ L q ) = 

n k kL 0 minjti ,ILq} 


= 5(0) - V £ £ 

k= 1 1=1, 


(lt\ 


j (It 2 {the integrand of Eq.(14)}. (18) 

l (fc-l)L 0 (l-l)L 0 

Assuming now that possible correlations between 5(£ 2 ) and bi,bj under the integral are 
small, S(t) itself varies very slowly within one correlation cell (see, however, below), and 
making use of statistical properties of random fields, (i) different transversal components 
within one cell are independent random variables, and (ii) magnetic fields in different cells 
do not correlate, we can average Eq.([L8]) thus obtaining a finite difference analogue: 

kL 0 


tl 


Sn = 5(0) - V £ < b{ k > S k 


dt\ 


dt 2 cos 25 (t 1 — tfi) = 


= 5(0) - 2 n 


k= 1 


,sin 2 5Lr 


{k—l)Lo {k—l)Lo 


5 2 


E < ^> S k , 

k=\ 


(19) 


where we retained possible slow space dependence of the r.m.s. magnetic field value < 
b 2 lk >= - <b\ >. Returning now to continuous version of Eq. (|T9p we get 


S{t) = --^L 0 

the solution of which has a form 


4 0 _ sin 2 5L 0 


3 r {5LqY 


l 

Jdt' < & 2 (t') > 5(0, 


Q(+\ Qfn\ ) 4 2 r sin 2 5Lq 

S(t) = S(0)exp{--n 


f < (riVj > dt’ 


( 20 ) 


( 21 ) 


For 5 -C L 0 1 and constant r.m.s. < b 2 > we obtain the simple (5-correlation result Eq. (|T5|) . 
The last form, however, allows to evaluate at what part of the convective zone the ASFC 
really takes place for a given profile of the r.m.s. < b 2 >. 
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3. Linked cluster expansion and dispersion 


Another important issue is the problem of temporal dependence of higher statistical 
moments of P aa . As P aa enter the Eq.([|) for the number of events one should be certain that 
the averaging procedure does not input large statistical errors, otherwise there will be no 
room for the solar neutrino puzzle itself. In order to apparently investigate generic statistical 


properties of the model we make an additional assumption and put in Eqs.(|13|), ([l4D S equal 
to zero. This is not critical at least for strong random magnetic fields as it follows that for 
small S instantaneous eigenvalues of Eq.(|l3l) 


E±(t) = ±\lP + n 2 b 2 ± (t) « 


( 22 ) 


i.e. for typical realization the vacuum oscillation parameter should not substantially influ¬ 
ence the solution provided also that 5 -C Lq 1 (see Eq. ( |21]) ). The integral equation (pA|) 
then takes the form 
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s(t) = s(o) - V / du / 


(23) 


especially convenient for perturbative solution in powers of magnetic held or, equally, in 
powers of coupling constant 4/i 2 . 

The n-th order is proportional to the product of 2 n transversal components b±(ti) inte¬ 
grated over time with descending upper limits of integration t > t\ > t 2 > ... > f 2 n _i. After 
averaging and some standard combinatorics (cf. 
providing the linked cluster expansion 


25j) the perturbation series exponentiates 


f 00 (—AiP)m. 

< S(t) >= 5(0) exp ( £ M 2m (t) 


(24) 


where because of isotropy we retained linked higher moments of even order only, which can 
be evaluated within our spacially piece-constant model of random magnetic fields as 


M2 mif) 




2m 5 


(25) 


where \ 2 m - some coefficients, 0 < \ 2 m < 1 , A 2 = 1 , dehned by magnetic held statis¬ 
tics within one cell. For Gaussian distribution all A 2m but A 2 disappear and Eq. fl24|) 
coincides with the 5-correlator result Eq.(15). Approximately this is also true when 


fl = (iL 0 ■ \J< b' 2 ± > 1 due to fast convergence of the sum in Eq. (|24|) . For chosen values of 

magnetic helds and correlation scale this parameter does not exceed ~ 0 . 2 , thus justifying 
the 6 - correlation estimates. It is also plausible (though we have not proved that correctly) 
that once we adopted the cell model with physically constrained from above random mag¬ 
netic helds, the resulting approach of < S(t) > to its asymptotic value should be always 
exponential leaving no room to such effects as intermittency. That is < S(t) > is simply 
decaying with time tending to zero. Here we illustrate such a behavior for a type of statistics 
differing from Gaussian. 


12 














Stochastic twist. Let transversal random magnetic fields in every correlation cell are 
constant by modulo, < b\ k >=< b\ >= const , differing from each other only by random 
direction of vector bj_ k . Then all even orders of the field entering Eqs.(|24]), fl25|) are constant 


too and \ 2 m — 1- The sum in Eq.fl24|) is easily performed and we obtain 


< S(t) >= 5(0) exp j — — sin 2 , 

l Ln ) 


(26) 


where = n\b±\L 0 . We note by passing that for large magnetic fields satisfying rather artifi¬ 
cial conditions = 7m, n — 1,2,..., there should exist an effect of resonance transparency, 
when the polarization vector 5 performs an integer or half-integer number of turnovers 
within one cell. Otherwise the behaviour is exponential again. 


To estimate the < S 2 (t) > behaviour we use the same procedure, write down the iterative 
solution for S(t), square it and then average. The final result again has an exponential form 
like ( |23|). As it turns out that fi < 1 we confine here only with the Gaussian statistics 
(valid in this case). Then, 

< Sit) >= 5(0)e _ri , < S 2 (t) >= (l + e -4r< ) , (27) 

where T is defined in Eq.([H]) and we also rewrite the expression (|15|) for < S(t) > . It follows 
from Eq. ([27]) that < 5 2 (0) >= 5 2 (0), and for t —> oo the exponent dies out and < S 2 (t) > 
tends to its asymptotic value < 5 2 (oo) >= |5 2 (0). For dispersion er| we then obtain 


= 7< S 2 (i) > - < S(t) >2 = h£) (! ^ e - 2 r< 


and correspondingly for <jp 


a P = P&t) > - < pjt) >2 = (1 _ e -2rt ) 


(28) 


2x'2 


(29) 


Taking into account that 


< f.,M >= < P ” (0) > (l + e- r ‘) , (30) 

we have that relative mean square deviation of P ee from its mean value tends with t —> oo 
to its maximum asymptotic value 


epjt) 

< Peeit) > 


71 ~ °' 707 ’ 


oo, 


(31) 


irrespectively of the initial value P ee (0). It is interesting that this asymptotic value is in a 
qualitative agreement with the result |26| , 


crpjt) 

< Pee[t) > 


~ 0.58, 


t —> oo, 


(32) 
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that was obtained by numerical simulation of the MSW-effect in a fluctuating matter density. 


We can now estimate the influence of the random magnetic field on the process of the 
neutrino conversion. Let suppose that the effective thickness of the part of the convective 
zone of the Sun carrying the large random magnetic field is about At ~ O.15P 0 and /i = 
10 ~ u /j,b- Then T sa 1.6 x 6/100 kG, s -1 . Hence, random magnetic field with the strength 
about 100 kG is strong, mixing parameter K = 1 — e~ rAt ~ 0.8. Random magnetic field 
b ~ 50 kG is medium, K sa 0.33 and 6 ~ 20 kG is weak, K sa 0.06. These estimations are 
the reasons to adopt the strength of the random magnetic field for the computer simulation, 
b = 50 kG and b = 100 kG. 


The estimation (pl|) of the r.m.s. deviation of P ee and other P aa is true, evidently, only 
for one neutrino ray. Averaging over N independent rays lowers the value ( |3T| ) in \/N times. 
That is for our case of N « 50 rays we get that maximum relative error should not exceed 
approximately 10%, thus justifying the validity of our approach. For smaller magnetic fields 
the situation is always better. 

To conclude this section, it is neccesary to repeat that the above estimate Eq. (|3l|) 
indicates possible danger when treating numerically the neutrino propagation in noisy media. 
Indeed, usually adopted one-dimensional (i.e. along one ray only) approximation for the 
(4 x 4) master equation Eq. (|ll|) or (2 x 2) Eq. (|T3|) can suffer from large dispersion 
errors and one should make certain precautions when averaging these equations over the 
random noise before numerical simulations. Otherwise, the resulting error might be even 
unpredictable. 


4- Interpretation of the influence of random magnetic field as a random walk over a circle 

Here we show that there exists a simple way to explain the behaviour of the mean value 
and dispersion of P ee (t). To apparently account the normalization P ee (t) + Pppit) = P ee (0), 
we can introduce a unit circle where a representative point parametrized by the angular 
variable 0(f) performs a motion, while P ee (t) = P ee (0) cos 2 0(f), Pp,p,(t) = P ee (0) sin 2 0(f). 
Then expressions for S and S 2 take form S(t) = P ee (t) — Pp.p,(t) = S'(O)cos20, S 2 (t) = 

S 2 ( o) 

—-— {1 + cos 40} . Suppose now that 0(f) is a realization of a Gaussian random process 
with the dispersion D</,(£) = fit and probability density 


P(4>) = 


y/2n/3t 


exp 


2 fit 


Simple calculation shows that 


< cos af> >= exp 


a 2 fit 


(33) 


From Eq.([33]) it follows that 


< S(t) > 
< S 2 {t ) > 


S{ fi)e~ 20t 

S\ 0 ) 




-m 
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‘-H U J /1 - 4 / 31 \ 


(34) 


Comparison of Eqs. ©, (|28|) with ( |34| ) shows that if we adopt 2f3 = T, these expressions 
become identical. Hence, one can treat the master equation ([i~3|) with the random magnetic 
held as a random walk over a circle with the dispersion of the angular variable proportional 
to the product of the path along the convective zone 5t, mean squared magnetic held < b'j_ > 
and correlation length L 0 , coefficient of proportionality depending on the neutrino magnetic 
moment squared, see Eq. (|T 6 |) . This provides a simple way to express the meaning of the 
Eq.(^ 0 |). At the starting moment 0(0) possess probability density concentrated at point 
0 = 0, and -P e e(0) has the definite value governed by the initial conditions. If t —> oo, the 
random walk over a circle provides an asymptotically uniform probability density, P^it = 
Too) = 1/27T and < P ee (t = Too) >= P ee (0)/2 because the average value of cos 2 0 is equal 
to 1 / 2 . A more detailed study of possible interrelation between neutrino conversions and the 
random walk will be reported elsewhere, here we only add that an explicit computation of 
the third and fourth moments of S conhrm the above identification. Our results also confirm 
the conjecture |5|] where with account of only one component of the transversal magnetic 
held (scalar case) it was shown that the resulting behaviour of < P ee > can be interpreted 
as a brownian motion of an auxiliary angular variable over a circle. 


D. Computer simulation of the master equation 

Substituting a random realization of magnetic helds along one ray (see curve 3 in Fig. |2] 
for b x (r) []) into the master equation Eq. (ED we fold the solution for four wave functions, 
v a (t) =| v a (t) | exp (ia a (t)), from which all dynamical probabibilities P aa obeying the con¬ 
dition Eq. (]]]) are derived. Then we have repeated such procedure with the solution of the 
Cauchy problem for other configurations of random magnetic helds (along other neutrino 
trajectories). 

After that, supposing that due to homogeneity the intensities of partial neutrino huxes 
are equal to each other, i.e. $f^(£, rf) = = const , we obtain the number event spectrum 

Eq. (Q) that depends on the product 

,n) = *(») x (P„>. 

Here is the integral flux of neutrinos of the kind ”i” |l^ and ( P aa ) are the mean 
probabilities that are shown in Figures [I] [J and in Figure |7| for particular cases of the 
maximum r.m.s. amplitudes yj ( b 2 ) = 50 kG and (b 2 ) = 100 kG correspondingly. 

For comparison of ASFC with the RSFP case we substituted regular magnetic helds 
Bo = 50 ,100 kG (with prohles shown in the same Fig. [| ) into Eq. (0) and calculated 
four probabilities P aa ■ The transition probabilities for antineutrinos P ss (appeared for large 
mixing angles in cascade conversions through the steps: (a) v e L —> within convective 


1 In general, for any random field realization the component b y (r ) has a different profile (along 
r) shown in Fig. [2] for b x (r) only. This difference in profiles was taken into account while equal 
amplitudes, \J(b x ) = f/(& 2 ) = \J(b 2 )/3 , were assumed too. 
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zone and (b) D^r —> v eR in solar wind via vacuum conversions) are shown for the cases 
B = 50 kG and B = 100 kG in the Figures |||9]. 

In Figures |3|-|9] along y-axis we plotted the dimensionless parameter 5 = 8.81 x 
10 8 Am 2 (eU 2 )/P(MeU). For instance, for boron neutrino energy E ~ 8.8 MeV and 
Am 2 = 10~ 8 eU 2 we find 5 ~ 1 . 

All four probabilities in the cases Figs. |3|-[7] obey the unitarity condition Eq. (H) for 
the same parameters sin 2 29, 5 = Am 2 /AE and for a given b rms . This can be viewed as a 
check of our numerical procedure because the unitarity condition was not apparently used 
in simulation. 

In Figs. [10]. [TT] we plot allowed parameter regions Am 2 , sin 2 29 found from reconcilement 
of all experimental data [lj, [|l 6 l for the ratio DATA/SSM ± lcr in the case of regular 
fields F>j_ = 50, 100 kG and in Figs. |T^, [L3] we present analogous results for random fields 
( b 2 ) = 50, 100 kG. In order to find these regions we have calculated number of events 
in Eq. (jTJ) using corresponding numerical solutions of the evolution equation Eq. ( |TT| ) with 
varying Am 2 , sin 2 29 devided by the number of events derived in SSM for P ee (E ) = 1. 


In Figs. |T0| - [T4| we also show Am 2 , 
antineutrinos v eR in SK |9fl. 


Sill' 


29- regions excluded from non-observation of 


IV. DISCUSSION OF RESULTS 


From Figs. m it apparently follows that there is a band in the parameter region 
Am 2 ~ 10~ 7 — 10 -6 eV 2 that separates MSW and magnetic held scenarios. In Figs. HH3 
we presented the allowed parameter regions to reconcile four solar neutrino experiments 
along with the SK bound on v eR . One can see that without the latter bound there exist two 
commonly adopted parameter regions, the small mixing and the large mixing ones, currently 
allowed (Figs. [L 2 ],[l 3 ]) or excluded (Fig. 0) in dependence of the values of the r.m.s. magnetic 
held parameters, see subsections IV.B and IV.C. But hrst of all we would like to discuss a 
” strange” dependence of the allowed from the SK experiment low mass difference regions 
( Am 2 < 1CU 8 eV 2 ) on the strength of the regular magnetic held. 


A. ” Paradox” of strong regular magnetic fields 


We find a surprising from the first sight result for the strongest field B = 100 kG: (i) the 
probability P gg decreases relative to the case B = 50 kG (compare Figs. 8,9); moreover, (ii) 
the well-known parameter region Am 2 ~ 10 ~ 8 eU 2 allowed here from SK data in the case 
of regular field B = 50 kG (shown in Fig. |K]) and recoinciled with other experiments in || 
vanishes in the case of more stronger field B = 100 kG (see Fig. pT|) . 

One can easily explain such peculiarity considering the analytic form of neutrino conver¬ 
sion probability in the convective zone (n eL —> v^ R ), 


(2jigx l 

(1/- Acos29) 2 + (2 /jB ± ) 2 


Sill 


\J(V - A cos 2 9) 2 + ( 2 /nBx) 2 -^-^ 


(35) 


This expression is valid for constant profiles and can be used to qualitatively estimate 
the effect treated here numerically for changing profiles. Here Ar = O.15P 0 is the 
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effective half-width of the convective zone where magnetic field is strong (see Fig. ^|); 
V = y/2 G F (p(r)/m p )(Y e — Y n ) < 1.6 x 10~ 14 (F e — Y n ) eV is the neutrino vector po¬ 
tential (V — V e — Vp) above the bottom of convective zone r > 0.7 R & (see Eq. ( 0 )); 
2[iB± ~ 5 x 10 15 B b0 eV is the magnetic field parameter for p = 10 11 ps and magnetic 
field strength normalized on 50 kG. 

The mass parameter A = 10 _ 11 Am 2 /2(E/MeK) eV is much less than parameters above 
for SK energies E > 6.5 MeV in the mass parameter region Am 2 < 10 ~ 8 eV 2 , or for 
A m\ < 10~ 3 . 

Since A m 2 {eV 2 ) = 4 x 10 -9 — 10 _1 ° is negligible and the resonance condition V = 
A cos 2 9, i.e., 


Am\cos26 

2E/MeV 


1.6(n - U) exp(-lO.54r/i? 0 ) , 


(36) 


is not fulfilled for corresponding low mass region, we conclude and check directly from the 
analytic formula Eq. (|35D that in the case B± = 50 kG the non-resonant large conversion 
probability P efi reaches a maximum due to accident coincidence of the phase in propagating 
factor with \jv 2 + (2pBj_) 2 Ar /2 ~ 7 t/ 2 . 

This is in contrast to another allowed region Am 2 ~ 10~' eV 2 shown in the same Fig. 
TO] where really the resonance Eq. (]36|) and RSFP take place. 

Notice that the oscillation depth (2pBj_) 2 /({V — A cos 2 6) 2 + (2pBj_) 2 ) reaches unity 
in Eq. (|35|) for extreme fields B± —> oo and the resonance V = A cos 26 (p 6 |) is almost 


irrelevant to the case B±_ = 100 kG. Similar behavior of RSFP was studied also in [[27]]. 

Moreover, for the case B± = 100 kG the phase in propagation factor reaches ~ 7 r for small 
mass parameters, Am] ~ 10 ^ 3 — 10 -5 resulting in zero (~ 0) of the conversion probability 
Eq. (fJH) while the survival one remains at the level P ee ~ 1. This is the reason why 
right-handed v e R are produced in a less amount for the case Rj_ = 100 kG than in the case 


B ± = 50 kG. 

This is also the reason why the low mass region Am] ~ 10 -3 — 10 ~ 5 which is allowed 
from SK data for B± = 50 kG (in Fig. (T0|) vanishes in the case of Bj_ = 100 kG (see Fig. 

ED- 

For low magnetic fields (we built but do not present plots for regular field B± = 20 kG) 
only resonant v e L —*■ iym-conversions take place in convective zone for Am 2 ~ 10~' eV 2 as 
well as MSW conversions for ’’large” mass parameters 10 -6 —10 ” 5 eV 2 . We can explain why 
this happens as follows. 

These fields are not too large to suppress a resonance in oscillation depth like for B ~ 
100 kG. On the other hand, for low magnetic fields RSFP occurs as strong non-adiabatic 
conversion for Am 2 ~ 10 -8 eV 2 in contrast to the adiabatic RSFP in the case Rj_ = 50 kG 
for the same mass parameter region. 

Therefore, this mass region is excluded from SK data in the case of low magnetic fields 
(for B = 20 kG) and the triangle region Am 2 > 10'' eV 2 which is common for all cases is 
allowed only. 

To resume the regular field case, we conclude that parameter regions allowed from SK 
data are close to analogous ones obtained in || for small mixing angle band, sin 2 26 < 0.25. 
However, in || exclusion of large mixing angles for small Am 2 relevant to RSFP comes 
from fit with other (radiochemical) experimental data while in our case the corresponding 
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region in Fig. 10 (where P g g occurs to be noticeable) is excluded from non-observation of 
antineutrinos v eR in SK and Kamiokande experiments || (see below). 


B. The SK experiment bound on electron antineutrinos and allowed parameter 

region 


In order to get antineutrino flux less than the background in SK we should claim |J 


$ S R SM (E > 8.3) til\ B (E)P s JE)a D (E)dE 4 2 . 

$ Pe (P > 8.3 MeV) = —— - K _ ’ y ’ -< 6 x 10 4 crrVV 1 


f s '* 3 \ B (E)a p (E)dE 


(37) 


where 

a p (E) = 9.2 x 10- 42 c7n 2 [(E - 1.3 MeV)/ 10 MeV} 2 

is the cross-section of the capture reaction u e p —> ne + ; $ B SM (E > 8.3) = & B J 3 5 3 X B (E)dE = 
1.7 x 10 6 cm~ 2 s is the SSM boron neutrino flux for chosen threshold of sensitivity E > 
E 0 = 8.3 MeV [|. 

Deviding the inequality Eq. 0 by this flux factor we find the bound on an averaged 
(over cross-section and spectrum) transition probability v eL —> V eR (via a cascade), 


$ Pe (P > 8.3 MeV) / 8 15 3 A B (E)P- eB (E)a p {E)dE 

& b sm {E > 8.3) / 8 15 3 X B (E)a D (E)dE 


or < 3.5% . 

We calculated boundaries in the parameter regions Am 2 , sin 2 26 where the inequality 
Eq. (^) is violated in order to exclude these regions from the allowed ones shown in Figs. 

EM- 

We do not find any violation of the bound Eq. (|38|) for low magnetic fields, B , b < 20 kG, 
both for regular and random ones. Thus, we conclude that in the limit B, b ^ 0 (for instance 
for the MSW case) the hole triangle region seen in Figs. 10-13 is allowed from the SK data. 

However, for strong magnetic fields such forbidden parameter regions appear in different 
areas over Am 2 and sin 2 29 for different kind of magnetic holds (regular and random ) . 
Moreover, for strong regular magnetic fields (B = 100 kG) such dangerous parameter regions 
vanish (see subsection IV.A above) while the stronger a random held would be the wider 
forbidden area arises. 

The bound Eq. (|38l) is not valid for low energy region below the threshold of the an- 
tineutrino capture by protons, E <1.8 MeV. 

Therefore we can expect in future BOREXINO experiment a large contribution of an¬ 
tineutrinos if neutrino has a large transition magnetic moment p ~ 10~ ll p B . This is due 
to re-scaling of the mass parameter Am 2 that is not hxed and decrease of neutrino energy 
that will be measured in BOREXINO. This prediction will not be in contradiction with the 
bound Eq.(|38D for hard neutrinos 

Really, the higher the neutrino energy (or the smaller the parameter 5 = Am 2 /4P) 
the less the probability {Pee) — {^tR^eR) occurs. One can see from Figs. |] r |] that for a 
wide region of the mixing angle sin 2 26 and, for instance, for the neutrino mass parameter 
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Am 2 ~ 1CU 8 eV 2 this probability is changing by 5-8 times: having maximum P ss ~ 0.25 — 0.4 
at low energies E < 1 MeV for sin 2 26 > 0.06 and the minimum (P g e 0.05) at the energy 
region E > 5 MeV. 


C. The MSW parameter region and correlation length of random fields 

For ” large” Am 2 ~ 4 x 10“' — 10 -5 eV 2 the MSW conversions without helicity change, 
^eL —>■ VfiL: take place under the bottom of the convective zone as one can see from the 
resonance condition Eq. (|3(|) where the neutron abundance Y n (neutral current contribution) 
should be omitted. 

Then left-handed neutrinos {y e L cross the convective zone with random magnetic 
fields and can be converted to right-handed components via the ASFC-mechanism. From 
Figs. 4, 7 it follows that the more intensive the r.m.s. held ^ ( b 2 ) in the convective zone 
the more effective spin-flavour conversions lead to the production of the right-handed P e R, 
D fl ji~ antineutrinos. This was also proved analytically in subsection III.C for small mixing 
angles both for h-correlated fields, Eq. ([T5|) , and for correlations of finite radius, Eq. 0). 

Only ’’large” mass region Am 2 > 3 x 10 6 eV 2 — 10 5 eV 2 survives as a pure MSW 
region without ASFC influencing the left-handed components since our correlation length 
Lq = 10 4 km chosen as a mesogranule size corresponds rather to finite correlation radius 
(see subsection III.C.2). Really, the dimensionless parameter x = 8L 0 is too large for these 
Am 2 and the ratio sin 2 x/x 2 ceases in Eq. (0)- Therefore, there are no ASFC for such mass 
parameters and Lq = 10 4 km. 

Variation of the correlation length L 0 . To check our approach we treated numerically a 
more stronger r.m.s. held b = 300 kG retaining invariant b 2 L 0 = const for a granule size 
Lq ~ 10 3 km, i.e. we considered small-scale h-correlated random helds. In this case a lot of 
v e R appear even for the typical small-mixing MSW region Am 2 ~ ICE 5 eV 2 excluding this 
and the whole triangle region at all from non-observation of v e R in SK. This immediately 
follows from analytic formula Eq. ( |2T| ) since for L 0 = 10 3 km limit (lim x ^ 0 sin x/x —> 1) is 
fulfilled there and the auxiliary function S ceases enhancing P Uf , L ^D^ R - 

The same situation remains even for a more realistic random magnetic held strength 
b = 100 kG with Lq ~ 10 3 km, see Fig. [T| where we reconciled all solar neutrino experiments 
with account of the bound Eq. (|38l)). 


V. CONCLUSIONS 


If antineutrinos v e R would be found with the positive signal in the Borexino experiment 
0 or, in other words, a small-mixing MSW solution to SNP fails, this will be a strong 
argument in favour of magnetic held scenario with ASFC in the presence of a large neutrino 
transition moment, /a ~ ICC 11 /!# for the same small mixing angle. 

There appears one additional parameter for the ASFC scenario comparing with the RSFP 
solution |0] to SNP. This is the correlation length of random magnetic helds L$ which varies 
within the interval Lq = 10 3 — 10 4 km in correspondence with typical inhomogeneity size 
(of granules and mesogranules) in the Sun. 
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The probabilitues P aa sharply depend on the correlation length Lq that might allow (in 
the case of v eR registered ) to study the structure of solar magnetic fields. 

Another regular magnetic field parameter /j,B± [0 changes to //\J(b 2 ) where the r.rn.s. 

magnetic field yj ( b 2 ) was treated here in the same interval b = 20 — 100 kG as for usual 
estimates of regular (toroidal) magnetic held in |7],§]. 

Our main assumption about more stronger random held is based on modern MHD models 
for solar magnetic helds where random helds are naturally much bigger than large-scale 
magnetic helds created and supported continuously from the small-scale random ones l|23l . |24|| 
(see subsection III.B). The ratio b = \JR m B |[23| . |24|| with large magnetic Reynolds number 
R m 1 means that in RSFP scenarios (including twist held model) the value of the regular 
large-scale held B± was rather overestimated. 

Thus, if neutrinos have a large transition magnetic moment p9[ their dynamics in the Sun 
is governed by random magnetic helds that , hrst, lead to aperiodic and rather non-resonant 
neutrino spin-flavor conversions, and second, inevitably lead to production of electron an¬ 
tineutrinos for low energy or large mass difference region. 

The search of bounds on // at the level /a ~ 10~ n /X-B in low energy v e-scattering, currently 
planning in laboratory experiments [Q, will be crucial for the model considered here. 

We would like to emphasize the importance of future low-energy neutrino experiments 
(BOREXINO, HELLAZ) which will be sensitive both to check the MSW scenario and the 
P e ij-pro duct ion through ASFP. As it was shown in a recent work [[□]] a different slope of 
energy spectrum prohles for different scenarios would be a crucial test in favour of the very 
mechanism providing the solution to SNP. 
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70000 km 


radiative convective zone 




FIG. 1. Geometry of neutrino trajectories in random magnetic fields. The finite radius of 
t'e^-neutrino source, the solar core radius O.I-Rq is shown. 


23 























z 


FIG. 2. Profiles of regular and random magnetic fields in the solar convective zone. One 
realization of random magnetic field component b x (r) (for one neutrino ray) is plotted. b y (r) 
component which has a different profile is not shown. 1 - regular magnetic field profile with 
B max = 100 kG, 2 - the same with B max = 50 kG, 3 - random magnetic field profile with 
brms = 50 kG. 
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FIG. 3. The averaged (over 50 rays) survival probability P ee = {v* eL v e i) in dependence on 
<5 = A m 2 /AE = 8.81 x 10 8 A m 2 (eV 2 )/E(MeV) and sin 2 29 - parameters for b r . m . s . = 50 kG and 
correlation length Lq = 10 4 km. 
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FIG. 7. As previous figure. The average probability Pgg = (i'eRi'eR) for = 100 kG and 

correlation length Lq = 10 4 km. 
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FIG. 8. The probability Pgg = (p*Rp e R) in dependence on 5 = Am 2 /AE and sin 2 20-paranreters 
for the regular field Bq = 50 kG. 
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FIG. 10. Am 2 (eV 2 ), sin 2 29 -plane for the ratio DATA/SSM in the regular field Bq = 50 kG. 
The isoline DATA/SSM = 0.368 ±0.026 allowed from the SK experiment is shown by the dashed 
line. The SK bound on v e R Eq. (|3^) is plotted too separating the excluded area. 
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FIG. 12. The reconcilement of all experiments. The isolines with dashed bands represent 
the DATA/SSM ratio equal correspondingly to 0.368 ± 0.026 (SK), 0.509 ± 0.059, 0.504 ± 0.085 
(GALLEX+SAGE), and 0.274±0.027 (Homestake). The random field b rms = 50 kG and correlation 
length L 0 = 10 4 km. 
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FIG. 13. As previous figure. The random field b rms = 100 kG and correlation length 
L 0 = 10 4 km. 
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FIG. 14. As previous figure. The random field b rms = 100 kG and correlation length 
Lq = 10 3 km. 
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